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The goal of this paper is to provide a new analysis of the classical dynamics of 

Bianchi type I, II and IX models by applying conventional Hamiltonian methods in the 

language of Ashtekhar variables. We show that Bianchi type II models can be seen as 

a perturbation of Bianchi I ones, and integrated. Bianchi IX models can be seen, in 

turn, as a perturbation of Bianchi lis, but here the integration algorithm breaks down. 

This is an "interesting failure", bringing light onto the chaotic nature of Bianchi type IX 

dynamics. As a by product of our analysis we filled some gaps in the literature, such us 

recovering the BKL map in this context. 

I. INTRODUCTION 



The goal of this paper is to provide a new analysis of the classical dynamics of Bianchi type I, II and IX 



On 

cr 

1 ' 

&Qmodels by applying conventional Hamiltonian methods in the language of Ashtekhar variables. We show 
> 

k> ^at Bianchi type II models can be seen as a perturbation of Bianchi I ones, and integrated. Bianchi IX 
models can be seen, in turn, as a perturbation of Bianchi lis, but here the integration algorithm breaks 
down. This is an " interesting failure" , bringing light onto the chaotic nature of Bianchi type IX dynamics. 

Bianchi models can be traced back to the nineteen century at least in the abstract, but their use 
in physical considerations had of course to wait for the invention of General Relativity (GR) and the 
development of cosmology. Bianchi I was first discussed by Kasner and further by Taub [BJ, Misner H 
and Lifshitz and Khalatnikov [Q.The empty Bianchi II model was solved first by Taub in a work where 
he discussed empty spaces with a three parameter group of motion || . Misner reformulated homogeneous 
cosmology using Hamiltonian methods ||, after the general hamiltonian formulation due to Arnowitt, 
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Deser and Misner (ADM) 0. This led Misner to introduce the fruitful concept of minisuperspace || 
where the cosmological evolution was equivalent to the motion of a particle in a potential. The cosmological 
implications of Bianchi IX were treated by Misner || and by Belinski, Lifshitz and Khalatnikov (BKL) 
|IIJ . We should note that when we write Bianchi IX, we restrict ourselves to the empty, diagonal case. |H| 
Bianchi IX is a homogeneous solution that became popular as a generic model of the approach to the 
singularity, even in inhomogeneous models Q. Although this is a controversial point [fT2"| |, the approach 
of treating the inhomogeneous case near the singularity as a perturbation (where the inhomogeneity were 
taking un account by a development in spatial derivatives ) from a zero order homogeneous solution 
(Bianchi IX) has been fruitfull ||13|| . In any case, the approach to the singularity as viewed by BKL shed 
new light on the dynamics of this complex solution since it can for be appproximated as a one dimensional 
map. During most of its evolution, Bianchi IX can be seen as a sequence of Kasner solutions, with the 
parameters of one Kasner "era" being mapped into those of the next; Barrow proved this map to be chaotic 

0- 

Actually, this was not the first time that chaos was proven in an approximation of Bianchi IX. Chitre 



T5fl proved that the Misner approximation of Bianchi IX was chaotic analitically using a know theorems 



about geodesic flows and the hyperbolic plane (see the article of Misner for a summary). Zardecki 



I7| first made numerical integrations of the Einstein Equations describing the Binchi IX evolution in the 



BKL form and obtained results that seemed to agree with these analysis. However, further investigations 



by Francisco and Matsas |1| and Rugh |19| showed that the initial conditions posed by Zardecki did 
not obey the scalar constraint. This could introduce negative energy densities, preventing collapse at 
the singularity and setting up chaotic motions. Furthermore, Francisco and Matsas computed the first 
Lyapunov exponent and showed that it tended toward zero. Analytical results from Burd et al. p(J and 
Hobill |21[ showed that the Lyapunov exponents had to be identically zero. Pullin and Rugh [^2| p3| 
explained this result by proving that the use of different time parametrizations could produce different 
Lyapunov exponents. The situation became confused, as usual chaos indicators seemed unable to capture 
the "irregular " character of (deterministic) Bianchi IX. By example a first analysis using The Painleve 
test seemed to indicate that it was integrable |24] but this claim was contested shortly thereafter [25 ^6 



using the perturbative Painleve test. Again, these works strongly indicated that the Mixmaster universe 
is indeed chaotic but we ought to remember that neither the Painleve test nor the perturbative one have 
reached theorem status, that is, at this time, they can only be taken as chaotic behavior indicators. The 
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latest original indication toward the chaotic character of the Mixmaster was given by Cornish and Levin 
|27| . They uncovered a (multi-)fractal repellor both in a two parameter map, the so-called Farey map ||28| , 
and in the full continuous dynamics. Fractality could be good chaos indicator in the context of GR since 
they cannot be undone via diffeomorphism. In the exact case however, the fractal structure is obtained 
numerically; as experience dictates, this is a tool to be use with caution with Bianchi IX. 

The main technical tool of our analysis is the appeal to Ashtekar's variables, introduced in the paper 



of Ashtekar of 1986 [p9 |. Ashtekars motivation was the quantization of GR, but he himself worked in the 
application of this new formalism in Bianchi cosmologies p0|. This was not the first time however that 



Ashtekar's formalism was applied to Bianchi cosmology |31 | . The introduction of the tetrad and first order 
formalisms was achieved in |3^| . Since 1986, the ideas of Ashtekar blossomed and ramificated in various 
directions as shown by the extensive bibliography on this subject [|33] The classical part is treated in detail 



by Romano. A short and pedagogical introduction tailored to our needs can be found in Giulini. 



Bianchi I and II were treated from a different point of view by Gonzalez and Tate. Relevant 

or in more 



treatments of the classical part can also be found in the article of Manojlovic and Mikovic. [37 



details and background in the book of Ashtekar |38[ . 

Our departure point is that Ashtekhar variables allow the simplification of the hamiltonian structure 
of Bianchi models to the point were they can be attacked by the methods of classical perturbation theory. 
These are almost as old as Newtonian dynamics itself. One of its uses was in celestial mechanics, where 
direct integration is not possible if more than two bodies are involved. The classic perturbation method 
based on Hamilton- Jacobi technique shows divergences, namelythe infamous small divisors problems 



11]. It was Poincare who first realized the complexity of the solution in phase space [ 40 1 . In this vein, a 



great breakthrought was the Kalmagorov-Arnold-Moser (KAM) Theorem [ [Hfl . 

By now, the onset of chaos in weakly perturbed Hamiltonian systems is well understood. In this case, 
the Hamiltonian splits into two parts: the integrable part H Q and the (small) perturbation AH. Usually, 
one makes a canonical transformation to action-ang le variables (/, 6) . Then H = H (i) + AH (f, 6) . 
Unperturbed motion is restrained to tori on phase space and is either periodic or quasiperiodic (integrability 
means that there are 2N constants of motion, there being 

2N degrees of freedom. In the Hamiltonian case, N constants of motion are sufficient to ensure 
integrability since the other N follow trivially from Hamilton equations. These N constants of motion are 
given here by /, which represents the radius of these tori). Expanding AH (f, in Fourier series 
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AH (l,0) =J2 AH (?) ex P 



n-e 



resonances will appear when N -0(1^) = where N are specific values for the n, u are the unperturbed 
frecuencies and I{q is some specific values of the actions I. Trying to get rid of theses primary resonances 
leads to secondary resonances. By example a canonical transformation near one resonance will lead to a 
Hamiltonian of the form 



i 1 

where U{ = dH jj , Ali = — Jjo and ip = N ■ 9+ (f>. This is formally the Hamiltonian of the nonlinear 



pendulum which its own new resonance term and its structure of elliptic and hyperbolic fix points. That is, 
at a smaller scale in phase space, a whole new resonant structure appears |42| KAM theorem ensures 
that under sufficiently small perturbations, almost all tori are preserved, but those near the resonances are 
destroyed. Another usual feature is the existence of fixed point; the saddle points are of particular interest 
leading to the appearance of stable and unstable manifolds and of the so-called homoclinic (heteroclinoc 
if more than one saddle point are involved) chaos [f4l |. 

One difficulty with the Mixmaster dynamics is that periodic or quasiperiodic trajectories (in minisu- 
perspace) do not exist and there is no fixed point. The absence of periodic or cuasiperiodic solutions is 
due to the monotonous increase (decrease) of the overall scale factor Q = In (abc) where a, b and c are 
the scale factors for the three axes. Indeed, it has been claimed that this absence would preclude chaotic 



behaviour for the Mixmaster [43|. But the monotonic increase (decrease) can be easily separated from the 
behaviour of the other significant variables and so one can recover quasiperiodic motions pTf . Solutions 
corresponding to a finite number of bounces before the trayectories make the perfect hit and go directly 
down one of the three channels are forbiden by the dynamics |I3 . 



All this paper can be seen as an exercise in classical mechanics in the somewhat unusual context of 
General Relativity. Our goal is to investigate Bianchi IX using Ashtekar's formulation. We will solve the 
integrable Bianchi I and Bianchi II models using Hamilton- Jacobi. We will describe the Kasner epochs in 
this formulation and latter apply simple canonical perturbation theory viewing Bianchi IX as a perturbation 
of Bainchi II. We will identify exactly where this approach breaks down, and analize the reasons for this 
"failure". 



II. BIANCHI MODELS AND ASHTEKAR'S VARIABLES 
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A. Homogeneous Cosmologies 



Among the cosmological model of General Relativity, some are particularly interesting for their math- 



ematical simplicity and their physical interest: the so-called Bianchi cosmological models [TjJ. Their 
simplicity resides in the fact that the spatial slice of these universe is homogeneous. The number of free- 
doms then reduces drastically and Einstein Equations become ordinary differential equations. The metric 
is 

ds 2 = -N 2 dt 2 + q u (t)x I ®X J ; I, J = 1,2,3 (1) 

Often this metric is written using the Misner's parametrization qu = exp (— 2Q(t)) (exp 2(5{t)) IJ where 
Tr{(3) = .This metric is (left) invariant under (spatial) transformations generated by a certain group 
of symmetries that characterize each specific Bianchi model. As usual, the Killing vectors £j are the 
infinitesimal generators of the isometries on these spaces . The left invariant vector fields Lj (x ! Lj = Sj ) 
verify then 

£ ij L I = [£, I ,L J ] = (2) 
This also means that the Killing vectors are right invariant vector fields |47fl. The left invariant vectors 



are related to the structure constants of the (Lie) group via 

[Lj,Lj] = C k jjL k (3) 

The C K j j are the structure constants of the Lie group that leaves the metric qu invariant. The 
spacetime spatial metric, that is, the metric written in a coordinate basis, is given by q^ u = quxixt- 
General Relativity admits a Hamiltonian formulation. Since we are interested in homogeneous cosmologies, 
we would like to write a simpler Hamiltonian, one which is homogeneous from the start. The Hamilton 
equations obtained from this simpler Hamiltonian will be the correct ones (that is, the same as the one 
obtains using the full Hamiltonian of General Relativity and demanding homogeneity afterward 
only when the structure constants verify 

C K u = e IJL S LK (4) 
were ejj# is the antisymmetric tensor and S IJ is a symmetric tensor. These models are called class A 



The simpler exemple is Bianchi I characterize by C K u = 0. The metric written in the Kasner form 



is: 



ds 2 = t 2pi (dx 1 ) + t 2p > (dx 2 ) + t 2 ^ (dx 3 ) (5) 
with pi + p 2 + Ps = 1 and p\+ p\+ p\ = ^- It is useful to use the following parametrization for the p's 



l + u 

P2 = ^TTT (7) 



u (u + 1) 

Ps = —. —r 8 

u 2 + u + 1 

where we assume the following ordering p\ < p 2 < P3 and 1 < w < 00. 

Another case of interest is Bianchi IX, C K u = e K u . Most of the evolution of Bianchi IX can be seen 
as a succession of Kasner epoch, that is the metric can be approximated to great accuracy by the Kasner 
metric. This stem from the fact that the potential consist of exponentially rising wall and are thus almost 
zero otherwise. Approximating these walls as vertical, the change from one Kasner solution to another one 
is given by the famous BKL u-map 

!u n - 1 if u n > 2 

(9) 
(u n - if u n < 2 

The u n > 2 case are called epoch and u n < 2 era. Chaotic behavior, if it happen, would be confined in the 

era change. This could be resumed by the following map 

u N+ i = u N - [u N ] (10) 

the Gauss map, relating one era with the next. Barrow proved this map to be chaotic in . 



B. Ashtekar formalism 

We will make a brief sketch of Ashtekar formalism. We will work mainly using two basis for spacetime, 
an orthonormal basis {e a } and a coordinate basis {d a } . Greek indices refer to coordinate (sometimes 
called spacetime in this context) bases and latin indices to frames bases (sometimes called internal indices) 
. When taken from the beginning of the alphabet (a, (3, a, b, ... ) their range is {0, 1, 2, 3} whereas for 
the middle of the alphabet (/i, z/, ... ) their range is {1, 2, 3} . 

The orthonormal basis and the metric are related through 

t] ab e a <8> e b = g a fidx a <g> dx 13 
Vabe a a e b p = g a/3 
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Thus taking determinant ont both sides gives 

e 2 = (det (e a a )f = det{g a0 ) = g 

viewing e a a as a square matrix of order n .We will fix the dimensionality of spacetime n = 4 . Let's 

fix an orthonormal base j ^ e± / 4 -* ej j so that ^ 4 ^ej_ is normal to the hypersurface S. The dual co-tetrad is 
j(4) e ± (4) e i j We have 

4 = iVe ± + N l ei 
at 

where N is the lapse function and N % e{ is the shift vector field. The relation between the spatial 
coordinate and frame basis is given by 

The point of departure is the action 

S = J{e)#x+F$[A]eZ4 (11) 



where +F^[+^A] = 2d [a ( +(4) ^] ab ) + 2 +{A) A [a a |c| +(4) A /3] cf> are the self-dual two-form. +{A) A a ab is 
a self-dual connection related to the ordinary (meaning neither self-dual or anti-self dual) connection by 
(see appendix 1) 

ab = \ [^A a ^ - \e°» cd ^A a ^ (12) 



+(^A a ab = i* ( +(4) A Q afe ) 



Y cd ^) Aa ^ (13) 



To recover the true degrees of freedom, as in conventional General Relativity, one makes an ADM 
decomposition ( J7|, see also [^]for an easier introduction; in this context, see f38|). The result is 

S = J {-ie* jk d t (^ jk E»)A * 

+KF« \ U E» E% + iN»i jk Fi\ u E»}d 3 xdt (14) 

where a surface term was dropped and is the covariant derivative constructed with A^ iJ which 
lives exclusively in the spatial hypersurface £ 

DpVt = d,V t + A^ % 



e ufe is the totally antysimmetric tensor density and the following notation was introduced 

^ s ^=dsW (15) 

N = —N = det(ef)N (16) 

~ vn 

E^ =_L£ ej 7 is the projection of the tetrad in the hypersurface E. The action is thus written explicitly 
in the form / [pq — H] whereas the Hamiltonian is a sum of constraints, since the conjugate momentum 
to the lapse, shift and the time-component of the 4-connection are absent. Our variables are jfc and 
II' 1 - fe , the self-dual part of — ie % j k Ef . The constraint's equations then read 

D, (f jk E?) = (17) 
F l ,M = (18) 
e iJ k F k = (19) 

To write the constraints as above, we used the fact that the self-dual Lorentz algebra is isomorphic to 
the Lie algebra of complexified SO (3) . The isomorphism can be carried out by taking an internal vector n a 
with n a n a — — 1. Then every self-dual internal 2-form f ab can be characterized completely by its "electric" 
part 2f a i>n a which is an internal vector orthogonal to n a . Hence 

Al = 2iAi b (20) 
K = 211^ (21) 
F a ^ = 2iF ± \ u (22) 



and we thus have 



I V - 2d [tx A v] 1 + ^ ik A {il j A v] k (23) 



To connect this formalism with Bianchi cosmologies let us introduce an orthonormal triad and its 
co-triad basis E\ 

dl 2 = 5 lJ E i E j (24) 
Our main interest in the introduction of this basis is to use them to define the Ashtekar's variables 

(£f,4) 
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£f = y/det {q^ESt (25) 

4 = r; - ii^ (26) 

The tilde denotes a tensor density. The are related to the Ricci rotation coefficients and the K l 
to the extrinsic curvature. Note the presence of i.Bj construction, Ashtekar's variables (at least in this 
formulation |[49|| ) are complex. We thus deal with a complex extension of General Relativity. In the end 
we will need to apply some "reality conditions " to ensure that the (spatial) metric constructed from the 
triad is real. Moreover one should also ask that the time evolution does not make the metric complex, so 
that the time derivative of the metric should also be real. Usually these variables are fields, that is they 
depend on x as well as on time t . In the homogeneous case, they will depend only on t . We can expand 
the triads on the invariant basis as follow 

E? = E?L*\x\ (27) 

4 = 44 (28) 

were |%| = det(x^) • Note that in this notation E\ = E^dij and E{ = Ejiq IJ that is capital latin 
indices are raised and lowered using the invariant metric qjj and the lower case latin indices using the 
orthonormal metric 5ij We have the following relation 

E Ii E i J = q IJ det(E!) (29) 

where det(E-) = det(g/j) . 

We will be particularly interested in the diagonal case: E{ ~ b\ and similarly for A\. We will denote 
the diagonal variables as Ej and Ai respectively. The Hamiltonian constraint H = reads 



(-^C J JK + lm A l jA™) EfE? 



E e, 3k {-AjC 1 JK + e J jk AjA K ) E 3 E K (30) 

J,K 



We have the following relations 



r / = r / = - eNLP5pi ^ {j} nlp ~ 1 C PNi^j 5R (32) 

- - x^f NLI — f-r -~n 

' ti E!\2 NLI A INL 



A I = A l I = ^-5 n = T I -iK I (33) 

e, - m 



The last relation can be easily inverted to give 

EiE 2 E 3 



Qd 



El 



(35) 



III. BIANCHI I 

As an introductory exercice and to explicitly relate this formalism to the usual Kasner solution, we will 
solve first the Bianchi I case. It is the simplest one, since all C K u = . Thus, the Hamiltonian is given 
by 

H = ( A 1 A 2 E 1 E 2 + A 1 A 3 E 1 E 3 + A 2 A 3 E 2 E 3 ) (36) 

To find the equations of motion, in this case we could integrate quite easily the Hamilton equations. 
Instead, as a preparation fot the more complex case we will find a convenient canonical transformation 
that makes the integration trivial. Consider the following generating function that implement a change 
from the old (phase space) coordinates (A, E^j to the new ones (j3, Coj : 

S(A, u) = ui In A x + u 2 In A 2 + u 3 In A 3 - Et (37) 

Thus 

dS uji , . 

Ei= dA* = A* (38) 

^ = ^ = ln^-|^t (39) 

and 

H = (iU 1 LU 2 + LUiLU 3 + U) 2 U) 3 ) (40) 

implying that uji are constants of motion and the f3 % are linear in the time parameter. Inverting we find 

Ei = ui exp (-fii) exp [-(u 2 + u 3 )t] (41) 

E 2 = u 2 exp (-p 2 ) exp [-(coi + u 3 )t] (42) 

E 3 = u 3 exp (-fa) exp [-(uji + 0J 2 )t\ (43) 
10 



Then 



<?22 



<?33 



UJ 2 UJ 3 

U) 2 
0J\0J 2 

^3 



exp ((3i - (3 2 - AO exp (-2cji) t 
exp (/3 2 -fh-Pa) ex P (~2^ 2 ) * 
exp (/? 3 - /3i - /3 2 ) exp (-2u 3 ) t 



(44) 
(45) 
(46) 



Up to now, our solution is complex. We have to apply the reality conditions. This means that both 
the metric and its time derivative must be real. This is easily seen to be obtained if 



(47) 



where Q is real and 



t = Reft] +ilm[t] 
= t + i Im[t] 



(48) 
(49) 



to a constant. We then obtain 



Qu 

<?22 
<?33 



2"3 



n t n 2 



exp 
exp 
exp 



2fiit 
2^ 2 t 
2fi 3 t 



(50) 
(51) 
(52) 



where we write t = Im[t] and the pref actor can be absorbed upon a rescaling of the axes. Also, General 
Relativity imposes 



H = 



then to recover General Relativity, the Q 's should obey the following constraint: 



fiifi 2 + fiifi 3 + n 2 n 3 = 



(53) 



(54) 



For convenience, to connect with the usual analysis in term of bounces and eras later on ||10|| , let us 
parametrize the Qi as follows 



^i = Po+P+ + V3p^ 
^2 = Po + P+ - V3p- 
tt 3 =p - 2p + 



(55) 
(56) 
(57) 
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Using the Hamiltonian constraint, we find 



Pl+pi=Po (58) 



Let's introduce the angle 9 : 



Let us define 



oi = 



p + = p cos 9 (59) 
P- — Po sin 9 (60) 



Qi 



Qi + q 2 + Q 



3 



2 ^ cos^ + ^V) (6i) 



3 V2 V 3 / 

Q 



0"2 



Qi + Q 2 + Q 



3 



2/1 cosfa-^Y) (62) 



3 V2 V 3 J 

Q 



o- 3 = 
~ 3 V2 

Note that <7i + <r 2 + (j 3 = 1 and of + of + erf = 1. Redefining time 



Qi + Q 2 + Q 3 
2/1 \ 

cos 9 (63) 



t= InT (64) 

Qx + Q 2 + Q 3 K J 



We obtain 



9ii = 


Q2Q3^ 2(Tl 
Ql 


<?22 = 


Q 2 


933 = 


QiQ2^ 2<J3 
Q 3 



(65) 
(66) 
(67) 

Which, upon rescaling of the coordinates x 1 — > rr 1 = Q^f 1 Q 2 Q 3 x 1 , etc gives the well known Kasner's 
line element 

ds 2 = T 2CT1 (dx 1 ) 2 + T 2CT2 (rfx 2 ) 2 + T 2 ^ (dx 3 ) 2 (68) 
It is useful to divide the domain of the angle 9 en six sectors 

(j-l)|<0<j| ; J €{1,2,3,4,5,6} (69) 
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Another parametrisation of the <7j is 



—u _ 1 + tt _ u(u + 1) 

u 2 + u + l u 2 + u + 1 u 2 + u + 1 



where 1 < w < oo and we write a to indicate the ordering <j\ < <5" 2 < 0*3 . It is straigtforward to write 
the u = u{9) . The result is 

u= l-a,(^,/3) 
cost' — 1/2 
1 + cos 9 
cos (9 — 2n/3) — 1/2 

where the j even (odd) formulae are generated each from the preceeding one using u{9) — > u{9) = 
u(0-2tt/3). 

IV. BIANCHI II 

This model is characterized by C 1 JK = S^jk • The Hamiltonian is thus given by 

H = ( A 1 A 2 E 1 E 2 + A 1 A 3 E 1 E 3 + A 2 A 3 E 2 E 3 ) - ^A 3 E,E 2 (73) 
We introduce explicitly the Newton constant G. When this is done we have f|8| 

4 = h ^ " O 

F* ^ = 2d [fM A v] ' + Ge' ]k A [fi j A v] k 
As in the Bianchi I case, we begin with a canonical transformation generated by 



F 2 = ln 

thus 



(GA 1 ) (g " 1 \GA 2 Y g ^\GA 3 )^ 3) ] (74) 



TT I \ (fa — fll — fa) /»r\ 

H = ( uj x <jJ 2 + ijJ X oj 3 + oj 2 uj 3 ) - u x u 2 exp — (75) 

Gr 

Since every one dimensional problem is integrable and since the (3 appears in a very special manner in 
the Hamiltonian, we propose a second canonical transformation, this time generated by 

G 2 = fixh + (3 2 1 2 + (/3 3 - A - (3 2 ) I 3 (76) 

We will then go from (8, Coj — > fa, i) . That will effectively convert the hamiltonian to an equivalent 
unidimensional one 
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H 



{hh 



I 2 



hh - (h + h) h + 4 



exp 



73 

G 



(77) 



Again, with an eye on future complications, we obtain the last canonical transformation to reduce this 
problem to a trivial one by applying the well know Hamilton- Jacob i technique. The generating function 



that will finally resolve our system, S |5C 



S = W{ l3 , P 1; P 2 , E) + 7l Px + 72 P 2 - Et 



(78) 



is a solution of the Hamilton- Jacobi equation 



that is 



/ dW \ OS 



E = P 1 P 2 



' dW s 



dW ( dW" 
Pl P 2 -{P 1 + P 2 )^L + [^L 

9l3 V ^73 , 

the other factors Piand P 2 are constants. The solution is 

G 



7s 

exp- 



G , n „ v f dx G f 
W = — (Pi + P 2 ) / 7 r e / 

2 v 2J J x + 1 2 J 



G f V cx 2 + bx 4 
x (x + 1) 



(79) 



(80) 
(81) 



(82) 



where e is ±1 



and we define 



Now 



c^(P 1 - P 2 ) 2 

b = -AE 

a = 4(PxP 2 - E) 



(83) 
(84) 
(85) 



x = exp 



73 

G 



(86) 



ap 

-G 



- 1 



e In 



2y/aJR(x) + bx + 2a 



^4ac — b 2 x 



(87) 
(88) 
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where we write R(x) = cx 2 + bx + a to abreviate the notation. Thus 

/ (t + 6 3 ) \ 2^E^R{x) + 2a + bx 
fsfexpl-v^e — - — 1= 1 (89) 

with 

if = V4ac - b 2 (90) 
Note that if posses units of E (same as units of I 2 ) . Writing 

u = exp (-yfa<t^*)\ (91) 



and solving for x, we obtain 

" (ifu - bf - 4ac 



AaKu 

x = 77^ — n2 : ( 92 ) 



We can now compute i 3 

aw 



(93) 

= -i e ^f M -^^ M+4 ty (94) 



•97; 

-1 

T^ v ~ (if « + 4 7 ) (K u + 4A) 

where we use the definitions 



7 = (2P 1 P 2 -E) + (Pi + P 2 ) v/PiP 2 - E (95) 



A = (2P 1 P 2 — E) — (Pi + P 2 ) JPxPs - P (96) 



The next step is to find 9i and 9 2 . 



ft = f (97) 

= 7l + Gln ^i±i2l + ^ lnll (98) 

(Ku + 4/3) >/a v 7 

similarly 

^ 72 + Gln (f^±M + ^ lnM (99) 
(if u — 4a) va 

were we make the choice e = — 1 (the e = — 1 case leads to an equivalent solution) and we used the 
following definitions 
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a 



= (P 1 - P 2 ) ^P,P 2 -E-E (100) 



P = (Pi- P 2 ) V p iP 2 -E + E (101) 
We are now able to compute the original Ashtekar's variables. Let us then compute the Ei. 



E, = ^ [ Pl + -Va) exp ^ j exp [~^^ ) (102) 
Wo , 1 r~\ f 6 A ( A(* + 3 ) N 



E, = W [* + ^) exp ) exp [~^^ ) (103) 
F = 1 (W- 167 2 ) 

x exp f- f h+k) ) exp f- f * + (t + ^ (104) 



g ; ; r v v g 

Writing down the metric is straightforward, since qn = ElE ^ 3 , and similarly the other components. 

Up to now, all our work was made with the asumption that the variables were complex. To connect 
with ordinary general relativity, we have to ask that qu and qu should be real. But the A 1 should be 
complex since P — iK l . Let us then choose 

P x = iVx , PiGRe + (105) 
P 2 = ip 2 , p 2 e Re+ (106) 
E,K e Re + (107) 

The hamiltonian constraint forces us to choose 

E = (108) 

thus 



V^ = 2^P 2 (109) 



a = JP,P 2 (P 1 - P 2 ) (110) 



7 = v^tv^+v^) 2 ( m ) 

To understand qualitatively the evolution let us focus on the essential features of the dynamics. Without 
the multiplicative constants, we have: 

ft=«p(-g) 

£i=exp(-^) (113) 



Ku V V G 
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where u = exp (— ^faG x t) . In the limit t — > — oo , we have 

Ei = exp (- P -§) 
£ 2 = exp (-^) 



E 3 = exp — JL — — t exp 



G 



Pi +P2 

G 



(115) 
(116) 

(117) 



In the limit t —>■ oo , we have 



-Ei = exp 
E 2 = exp 



P2t 
G 
Pit 
G 

E 3 = exp f ^fi j ex P 



Pi +P2 



(118) 
(119) 

(120) 



These asymptotic forms for i£j have the typical aspect of a Bianchi I solution 

E 1 = exp[-(0 2 + 3 )t] 

E 2 = exp[-(n 1 + n 3 )t] 

£ 3 = exp[-(Oi + 2 )t] 

and 



(121) 
(122) 
(123) 



Ei = exp 
E 2 = exp 
£3 = exp 



-(Of + Of )*' 
-(Of + Of )t 



(124) 
(125) 
(126) 



thus we can approximate the evolution as a transition from one Bianchi I-like state to another (the 
famous "bounce" in the usual minisuperspace version of homogeneous cosmology, where the evolution is 
seen as the motion of a particle in a potential [H]]). The transition from one Bianchi I history to another 
is better understood if one redefines new 0,- and O^ such that 



Oi < 2 < O s 



of < o 2 v < o 



TV 



(127) 
(128) 



Then 
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af = ft 2 + 2fi 3 



-a 



(129) 
(130) 
(131) 



N 



Introducing a new angle 9 

pl = p» cos 9 n 
p^=p$sm9 N 

we can easily find the transformation law for the parameter 9 

N 4 - 5 cos 9 
cos 9 - 

sm9 N 



(132) 
(133) 



5 - 4 cos 9 
3 sin# 

5 - 4 cos 9 



(134) 
(135) 



This is of course essentially the BKL ]TI| analysis seen in a somewhat unfamiliar context. 



V. BIANCHI IX 

In this case, we will see that a straightforward application of classical perturbation theory fails to find 
a canonical transformation that makes the Hamiltonian trivial and thus the integration possible. This 



failure is of course perfectly understandable in the light of recent developments concerning Bianchi IX ^5 
[ 26] |27]].The Hamiltonian is 



H = ( A 1 A 2 E 1 E 2 + A.AsE.Es + A 2 A 3 E 2 E 3 ) - 
1 



G 



{A z ErE 2 + A 1 E 2 E 3 + A 2 E Z E 1 ) 



(136) 



By inspection, one recognizes part of this Hamiltonian as the Bianchi II that was shown previously to 
be integrable. Using the same canonical transformation as before, namely equations ( |74]) and (0), we find 



H=(I 1 I, 



cxp 



73 

G 



(T j s T - (27i + 7a) (T j , T - (272 + 7s) 
(ii - J 3 ) J 3 exp (I 2 - J 3 ) J 3 exp - 



G 



G 



(137) 



With a further transformation 



S 2 = H/( 73 , Pi, P 2 , E) + lx P x + 72P2 



(138) 
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using ([lOg, [TOT], H and ||) one finds 



P2 + ^y/a 



'Ku + 4(3" 
Ku + 47 / 
fT-u — 4a ' 
Km + 47 



(139) 
(140) 



H = E 



1 (Km -4 7 ) f / 1 >v .„ , . -201 
^=^-^ I (p, + - VS) (*« - 4«) exp _ exp 



-2P,0 



2 U 3 



G 



(141) 



( P 2 + - v^J (Ku + 4(3) exp exp — -±± 



We know that General Relativity will require H = 0. Since we want to consider a weak perturbation, 
we can set E = in the perturbed Hamiltonian. 



o 



7 



P X P 2 (P - P 2 



ppf\/p +V^) 2 



w — exp I -2a/PP 2 



K = 4JPP 2 (P x - P 2 



Thus 



H = E + 



1 



P tM + JP 2 (P x -P 2 )u- 2P 1 JP 



+ \/Pi ( \/P 
+ 



P> u 



-1 



-201 -2P 2 fl 3 
exp — — — exp 



G G 
P~2 (JPi + \fK) (Pi -P 2 )u- 2P 2 (Jp'+^PtJ 



P 2 [JPi + JP 2 ) M" 1 



exp ■ 



-29 2 
G 



exp 



G 



We then have the following Hamiltonian 



H(P,9) = H (P) + AH(P,9) 



We are looking for a generating function F: 



F = 9- J + X(9,J) 



(142) 
(143) 
(144) 
(145) 

(146) 



(147) 



(148) 



(149) 
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so that 



H(P,9)^H (J) (150) 

up to first order. Thus F generates a canonical transformation from (O, to Jj .We want to 
integrate the system up to first order in the perturbation. We know that 

dF 2 dX , s 

p > = w = J > + w (151 > 
*=£H + f 



By Taylor 



H(P, 6) = H(J, $) + — (61* - <H) + — (P -Ji) + (e 2 ) (153) 
= H (J) + AH (J, $) + + O ( e 2 ) (154) 



Our case is even simpler since 



Thus 



If X is so that 



then 



that is 



d9 3 

which is the same as asking that 

dx 

or 



H = E(=P 3 ) (155) 



BX 

H(P, 6) = Ho(J) + AH (J, 0) + — + O (e 2 ) (156) 



H(P,6)^H (T) + O(e 2 ) (157) 



AH(fJ) + — = 0(e 2 ) (158) 



^- = -AH(fJ) + 0(e 2 ) (159) 



- A H( J, 6) + O (e 2 ) (160) 
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X = - J AH (J, 6)d6 3 + O 



(161) 



Thus 



X 



G 



-== exp 



G 



-2P 1 (yp[+vptf - 2 p 2 g 3 



cxp • 



G 



/ P~l(VP~l + VPt) 3 -2yfPi(yfP 2 -y/Pl)0 Z 

■ exp 



G 



exp • 



-2gi 
~G~ 



(162) 



+ 



ir (Pi _p 2) -2^(^+^)03 2P 2 (v^+v^) 2 _2P^ 3 
— — exp — — — exp ' 



'Pi 
(v'P~ 



"P~i(VP~i 



G 



-2y/P[(y/I\-y/F£)e 3 

exp 



G 



exp ■ 



-Pi 

-2£ 2 
G 



G 



is the generating function that would permit the integration of the whole system albeit to first order. 
But the generating function X diverges and is thus useless near the following planes 



P 1 = 
P 2 = 
Pi = P 2 



(163) 
(164) 
(165) 



The physical significance of these values is explained in the appendix 2. 

Let us analyse the P 2 = resonance. It is clear that the problem with X arises because the exponents 
lead to inverse powers of P 2 in the generating function; positive powers of P 2 in the pre exponential factors 
only improve integrability; thus, to analyze whether the system may be integrated or not, we may drop all 
powers of P 2 from the prefactors and write H as 

1 



H = E + -P? exp 
4 



E+-P?exp 



-29! 




-2P 2 6 3 


G 


exp 


G 


-26t 




-2P 2 e 3 


G 


exp 


G 



— 2\/ P 2 \/ P\Q 3 2\/P 2 \jP\Q 3 
-2 + exp — h exp - 



G 

' , 2^Y 2 ^P[B 3 

cosh — 1 



G 



(166) 



The typical (perturbative) route to chaos is seen by seeking second order resonances. Let's then 
integrate the more resonant term using the following generating function 

e6 3 



Get exp ( ^ j + s6 2 + K9 3 - f ex P 



-2- 



G 



- 1 



(167) 
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which allows us to make the change from (9 1 , 9 2 , 9 3 , P 1 , P 2 , E) — > (fa, fa, fa, £, ^0 • The new Hamil- 
tonian reads 



1 2 -2e^ 3 , 2Vefa 
H = K+-a 1 exp - cosh j - ^ 



a a 2 / -2e^ 3 N 
G^ 1+ 2F eXP ^ \ 



168) 



where 



r ^ n , 1 D 2 —1Q\ -2P 2 ^3 
2 1 6XP "g" 6XP — G — 



(169) 



While we seem to have isolated the resonance at e — 0, the new Hamiltonian has other secundary 
resonances. To bring them forth, let us attempt to integrate the new Hamiltonian to first orden using 



where 



Consider 



where 



Thus 



X 



AKd^ 3 
Ga 2 



Ae 



J exp (—z) cosh ^z\J A + B exp (—z) j dz 



A 
B 



a a 



a (fa 
e \ G 



a 
2e 



a 



2e 2 



J exp (—z) exp (^zA\Jl + B exp (—z)j dz 
/exj 



-1 + Ay/l + 5 exp (-z)J -2 
y exp (—z) exp ^— 2:^1^ -KB exp (—z)j dz 
J exp - (l + Ay/l + 5exp(-^)^) z 



1 



B = — = 



B 

A 



a (fa 
e\G 

a 



a 
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2, ($ - 



2c 



(170) 
(171) 



(172) 
(173) 



(174) 



(175) 



(176) 
(177) 



(178) 
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Let's analyse the divergence of I± . Let's tackle the case of J + first. 



k=0 



(179) 



where 



1 + A\Jl + Bexp (-z) ) z 

— [ x k exp (— x) (1 + B exp (— x)) k ^ 2 dx 
k\ Jo 

(l + 5exp(-A;)) fc/2 



180) 
181) 



The series diverges when A — 1 . We argue that A = 1 is a (simple) pole for this series, at least for 
some range of the initial parameters. To show this, we will construct an analytic continuation for the series 
(see appendix 3). Then 



k=0 



k=0 



A 



where 

l (VT+B) k -l 

FJB) < -B K - > 

We argue that A = 1 is a simple pole since the series 

OO , A \ k 

k=0 



are analytic in A = 1. Using equation |176| we encounter "resonance" for the values 



a (ipi a 



That is 



1 



1 



e=-P 1 ±-^P 1 2 -2a 



where we used the generating function 167 to obtain the relation between a,i[>i and Pi 



a 
G 



ih = Pi 



182) 



183) 



184) 



185) 



186) 



(187) 
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VI. CONCLUSION 



Using Ashtekar formalism, we treat the Bianchi IX Hamiltonian using the tools of classical perturbation 
theory. Namely, we treat Bianchi IX as a perturbation of the integrable Bianchi II. The failure to integrate 
Bianchi IX perturbatively did not come as a surprise in view of the numerous analytical and numerical 
evidences of its chaotic character. Our main purpose was to show how the use of Ashtekar variable simplify 
the perturbative analysis since the Hamiltonian (often refered to as the scalar constraint in this context) is 
notably simpler then in the usual one. As a by product of our analysis we filled some gaps in the literature, 
such us recovering the BKL map in this context. The reality conditions, which are often difficult to handle 
are easily deal with here using the known Bianchi II solution and the relation between the two formalisms. 

From the point of view of the larger framework of chaos and general relativity, the important point is 
to analyze the reasons for the "failure" of our naive approach to Bianchi IX. Our strategy was simply to 
apply to Bianchi IX, once rendered manageable by the translation into Ashtekhars variables, the most direct 
known ways to handle a Hamiltonian system. We treaded at first on known grounds, and not surprisingly 
found a quick success (providing a more complete solution to the Bianchi II case than previously reported, 
a tribute to the power of Hamilton - Jacobi methods). Thus encouraged, we faced the Bianchi IX problem. 

The expectation in facing a complex Hamiltonian dynamical system, but that can be divided into an 
integrable part and a "perturbation", is that either no resonances will appear, and then the system will 
be integrable by virtue of KAM theory, or else there will be resonances. If these are isolated, however, 
the dynamics is not yet chaotic, but rather the resonances appear as boundaries of regions of integrability. 
Chaos arises when resonances appear in layers, an infinite set of "secondary" resonances accumulating 
towards the primary ones. 

Because of the complex nature of Ashtekhar variables, it could not be expected that this analysis would 
translate in any straightforwrad way to our problem; rather, our aim has been to discern whether a similar 
structure to the familiar weak chaos appeared in our problem. Thus, instead of primary resonances, we 
found that naive perturbation theory (seeing Bianchi IX as a perturbation of Bianchi II) breaks down 
in channel runs, and also near the singular solutions at the end of the channels. The issue then became 
whether these were isolated singularities, or rather there were other breakdown points arbitrarily close to 
them. 

In order to find an answer, we zeroed on the P 2 = singular metric, by isolating the terms responsible 
for the breakdown of naive perturbation theory, and proceeded to the exact integration of this most 
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singular part of the Hamiltonian. If (to our surprise) no further obstructions to integrability had appeared, 
this would have been an indicator of a non chaotic nature of the whole Bianchi IX system; alas, our 
expectations held on, and we found that, arbitrarily close to the original singular solution, new singularities 
of perturbation theory appeared, as given by Eqs. (186) and (187). 

To analyze these conditions, we must remember that, from the point of view of the dynamics generated 
by the Hamiltonian K, both a and ipx are constants of motion. Thus what we found is, for any given 
value of the original momentum P 1; an hyperbola of trajectories where perturbation theory breaks down, 
approaching asymptotically the known singularity at e = when a — > 0. The P<i = e = singular point 
cannot be isolated, therefore, and we must expect complex behavior in a neigborhood of the corners of 
Misners triangular potential well. 

We see that the failure of our naive approach in fact is leading us directly to the sources of complexity 
in the dynamics. This bolsters the conclusion that Ashtekar variables are a useful tool in studying chaotic 
behavior in General Relativity. 

We are indebted to Luca Bombelli, Viqar Husain, Ted Jacobson, Jorge Pullin and Michael Ryan for 
enlightening discussions. This work was partially supported by Universidad de Buenos Aires, CONICET, 
Fundacion Antorchas, and by the Comission of the European Communities under contract Nr. C11*-CJ94- 
0004. 

VII. APPENDIX 1 



We will resume here some basic notions of differential geometry. Our notation follows [|lq| . First, given 
two n-forms A and a at some point in M : 

A = ^A ai ... an e ai A...Ae°" (188) 

71! 

o = V»i-«* eai A - A ^ ( 189 ) 
we defined their inner product, induced by g , via 

< A, a > = ^A ai ... an g a ^...g a ^a bl ... bn = \ ai ... an (r ai "' an (190) 

The following proposition follows readily 

\A* a =< \,a > e (191) 
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where e is the volume form on M induced by the metric 



e = y/gdx° A ... A c£r (n ~ 1 - ) 



e ai ...a n e ai A ... A e u 



nl 



(192) 
(193) 



where e ai ... Qn = n\5° [ai ...5^ n] l) . 

Consider the following string of equalities 



S 



J F ab A* (e a A e b ) 

\ J Rabcd (e c A e d ) A* (e a A e b ) 
1 f 

- J Rabcd < e c A e d , e a A e b > e 



Re 



(194) 



This show how to write the Einstein-Hilbert action in terms of the curvature 2-form F a b and the 
co-tetrads e a . The curvature F ab is defined as 



F a b = -Rabcd e c A e d 



and satisfy the second Cartan 's structure equation 



(195) 



du; a b + u a c Au c b = F ab 



(196) 



The u a b are the connection 1-forms, defined as 



,,, a — r a p 

U b — 1 bc e 



(197) 



The component of the u a b in coordinate base are called the Ricci rotation coefficients [47 



T a bc e c a dx a = u a \dx a 



(198) 



That is 



c -pa 
e a L be 



(199) 



It is straightforward to prove that uo ab = —tOb a [0. The following definition is useful 
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V a V = V a (V a e a ) 

= d a (K) e a + V a e b a V b e a 
= (d a (V a )-V c u a \)e a 
= {d a (K) + V c cu aa c ) e a = D a (K) e a 

Finally the action can be rewritten as follows 

S = J F ab A* (e a A e b ) 

= \J R abcd (e c Ae d ) A* (e a A e b ) 

= \j R aba pe a c e P d (e c Ae d ) A* (e a A e b ) 

= \J Raba^e 13 , < e c A e d , e a A e b > e 

= \j Raba^el (r,^ - 7fV C ) 6 

= / R Cd a ,e a A(e)d'x 
where d 4 x is a short notation for dx° A ... A dx^ n ~^ . 



VIII. APPENDIX 2 

The purpose of this appendix is to give the physical meaning of the resonances 

P 1 = 
P 2 = 

First let us write the metric: 

-1 (Kuf-l^ 2 ( 6 2 \ ( e 3 \ 



-1 (i^) 2 -16 7 2 / M ( 9p 
922 = IfifbG K~u 6XP \ ~G ) 6XP (" 2P2 



0_3 

G 



where 
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u = exp (-y/p[p 2 ^j (208) 

7 = v^^V^+V^) 2 ( 209 ) 

We can rewrite this in a more convenient manner. To recover real general relativity, we choose Pi and 
P 2 imaginary and write separately the real and imaginary part of #3. We finally find 



qn = " yW ^ 2G V ^ 7 ^ cosh (V^jjf J exp ^-2g J exp J (210) 

922 . z^|±^ cosh (^|) exp (_ 2 |) exp (_ 2P2 g (211) 

1 iVTTft 

* 3 " G^(7W^) (212) 

where we write 6> 3 instead of lm[0 3 ] to simplify the notation. If P 1 = P 2 , qn = q 22 . 
Recall that, in the standard Misner notation, 

qn = exp 2 (-ft + ft + v^ft) (213) 

g 22 = exp 2 (-ft + ft - v^ft) (214) 

Thus 

9n = <? 22 ft = (215) 

Viewing the evolution of the universe as a scattering process in a triangular potential well, this gives 
us a particle-universe running right into the right hand channel. 
Let's see now the case P 2 — > . In this case 

qn — > 00 (216) 

while the other two components of the metric go to zero. This correspond to a particle going into one 
of the left hand channel. Similar behavior occurs with P 1 — > 0. 

IX. APPENDIX 3 

Let's go back to our exact integral for I k : 

4 = 1 r x k exp{-x) (1 + B exp {-x)) h/2 dx (217) 
k\ Jo 
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We can write 



(1 + B exp (-x)) k/2 = 1 + Ul + Bexp(-x) - lj £ + Bexp (-x)) ' 

J i=o ^ 7 

If \B\ < 1 , then 

yl + 5 exp (-rr) < 1 + -5 exp (-x) 



(218) 



(219) 



and 



k-l 



k-l 



J2Ul + Be W (-x)) <J2{VTTB) 1 

7 i=0 

(v / TT5) fe - 1 



vTTb - 1 



(220) 
(221) 



Thus 



k= k\\l xk ex P (~ x ) dx + J x h exp(-x) yl + 5 exp (-x) - 1 ^1 + 5 exp (-rr)^ dxj 

1 / "°° *■ ' N l ^Jl + Bexp(-x) - 1 ^ ^l + Sexp(-x)y dx (222) 



= I f / x fc exp (-x) 



now 



roo 

/ x fc exp(— x) 
Jo 



k-l 



^/ / V 1 (VTTb) -1 l ... 

^l + J Bexp(-x)-l ]T ^l + Eexp(-:r)W:r<-E V / __ / — - / y k ex P (-y)dy 



i=0 



4 v 7 !^^ - 1 2 k Jo 

i (v / TTs) fc -i i 

4 y/T+B - 1 2 k 



Therefore 



where 



h = l + F k (B) 



i (VTTb) - 1 i 

FAB) < -B K t J ^ 

V ; ~ 4 2 fc 



(223) 



(224) 
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